Energy gaps, magnetism, and electric field effects in bilayer graphene nanoribbons 
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Using a first principles density functional electronic structure method, we study the energy gaps 
and magnetism in bilayer graphene nanoribbons as a function of the ribbon width and the strength 
of an external electric field between the layers. We assume AB (Bernal) stacking and consider both 
armchair and zigzag edges and two edge alignments distinguished by different ways of shifting the 
top layer with respect to the other. Armchair ribbons exhibit three classes of bilayer gaps which 
decrease with increasing ribbon width. An external electric field between the layers increases the gap 
in narrow ribbons and decreases the gap for wide ribbons, a property which can be understood semi- 
analytically using a 7r-band tight-binding model and perturbation theory. The magnetic properties 
of zigzag edge ribbons are different for the two different edge alignments, and not robust for all 
exchange-correlation approximations considered. Bilayer ribbon gaps are sensitive to the presence 
or absence of magnetism. 

PACS numbers: 71.15.Mb, 81.05.Uw, 75.75.+a 



I. INTRODUCTION 

Steady experimental and theoretical progress^ in un- 
derstanding the physics of single and multilayer graphene 
sheets, with and without substrates^ 2 - has attracted at- 
tention from the technical community interested in ex- 
ploring the potential^, of this truly two-dimensional ma- 
terial in electronics. Because graphene is atomically 
thin, it automatically scales channel-thicknesses to at- 
tractive values. Isolated single-layer graphene sheets are 
zero-gap semiconductors when extremely weak intrinsic 
and Rashba spin-orbit interactions (SOIs)^) are ignored, 
and have a room-temperature carrier mobility that is 
weakly 5,6 carrier-density dependent and higher than in 
known compound semiconductors. Because of the rela- 
tively weak SOIs, and also because the zigzag edges of 
graphene sheets have been predicted by theory 7,8 to be 
magnetic, there is also interest in graphene as a poten- 
tial material for spintronics. Initial experimental efforts 
in this direction have focused on injecting spin-polarized 
carriers from magnetic metals A 

Before graphene can be used as a replacement for or 
a supplement to silicon in CMOS circuits, it will be 
necessary to prepare graphene channels to have a high- 
resistance or the off state. The most obvious way to 
achieve off states similar to those of silicon is to in- 
duce an energy gap in the density-of-states of a graphene 
system. One possibility is to open gaps by tailoring 
sheet-substrate interactions. Indeed angle-resolved ph- 
toelectron spectroscopy (ARPES) spectra of single-layer 
graphene sheets on SiC substrates^ hint at this possi- 
bility and have been interpreted as showing gaps ~ 0.26 
eV. A more obvious route is to induce sizable quantum- 
size gaps by using narrow ribbon&Ji. One early exper- 
imental report of energy band engineering in graphene 
nanoribbons 3,2 has already appeared in the literature. A 
second strategj*i£ for opening gaps in graphene systems 
is to use a bilayer geometry and apply an external elec- 



tric field directed between the layers. Recent Shubnikov- 
de Hass cyclotron-mass measurements^ and transport 
studies^ 5 - in bilayer ribbons suggest that such an elec- 
tric field does indeed open up a gap. The transport 
measurements suggest that a gap of the order of meV 
opens up when an electric field of about 0.167 V/nm 
(which corresponds to applying 50 Volts across a 300 nm 
Si02 substrate) is applied across the bilayer. A distinct 
advantage of bilayer ribbons in nanoelectronics was re- 
cently suggested by the IBM group^. Their experiments 
hint at suppression of electrical noise in bilayer graphene 
channels comapared to the noise present in their single- 
layer counterparts, thus improving the signal-to-noise ra- 
tio. Recently, using a tight binding (TB) modeli 7 -, two 
families of zero-energy edge states were found in bilayer 
zigzag ribbons. 
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FIG. 1: (Color online) Schematic illustration of the two edge 
alignments we consider in bilayer graphene. (a) Q-alignment 
and (b) /3-alignment. The /3-alignment can be obtained from 
the a-alignment by shifting the top layer. 

In this article we use density-functional theory^ 
(DFT) based first principles methods 19 and tight-binding 
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(TB) models to explore the physics of energy gaps in 
graphene systems when the narrow ribbon and external 
field approaches are combined. We assume AB (Bernal) 
stacking and consider two types of edge alignments (Fig. 
[1]) for both armchair and zigzag ribbons, denoted as a 
and (3 alignments. (The top layer in /3-alignment is 
shifted with respect to that in the a- alignment.) We 
report on the ribbon-width and electric field dependence 
of bilayer gaps for both edge types (armchair and zigzag) 
and both edge alignments. Because zigzag edges in par- 
ticular have a tendency toward magnetism which, when 
present, can have a large impact on energy gaps, an im- 
portant element of this study is an assessment of mag- 
netic tendencies in bilayer ribbons. To the best of our 
knowledge there are no previous systematic theoretical 
studies of ribbon-width, magnetism and electric field ef- 
fects on the energy gaps in bilayer nanoribbons with a- 
and /3-alignments. 

The physics of graphene nanoribbon systems is en- 
riched by the theoretical possibility of broken symmetry 
states. The ferromagnetism predicted for single layer rib- 
bons with zigzag edges^^, which is related to 7r-electron 
edge orbitals, is one possibility. The edge states were seen 
in non-magnetic STM experiments^ -. Interestingly half- 
metallic magnetism has recently been predicted^ when 
an electric field is applied across a ribbon with zigzag 
edges and half-metallicity is argued to be sensitive to 
the nature of exchange-correlation potential^ 2 .. One first- 
principles DFT study of edge ferromagnetism in bilayer 
zigzag ribbons^ 3 -, which employs the local-density ap- 
proximation (LDA)2i for the exchange correlation poten- 
tial, predicts that a non-magnetic ground state appears 
when two single-layer graphene systems are stacked. Our 
DFT calculations also indicate that edge magnetism is 
less robust in bilayers than in single-layer systems. We 
have found that DFT predictions for the magnetic state 
of bilayer zigzag ribbon systems, in both edge align- 
ments, are sensitive to the particular semi-local approx- 
imation that is employed; magnetism appears when the 
generalized gradient approximations (GGA), PW91— or 
PBE962& are used, but not in the LDA. Even in single- 
layer graphene, one recent study2£ has indicated that 
edge disorder strongly influences edge magnetism and 
hence ribbon gaps. The possibility of a novel broken 
symmetry associated with non-local exchange effect, re- 
cently predicted using Hartree-Fock theory 2 ^, peculiar to 
Bernal stacked graphene which would influence ribbon 
gaps, is not considered here since the physics behind it 
would not be picked up by any general-purpose exchange- 
correlation approximation. For all these reasons, the ex- 
perimental predictions from the present DFT study of 
bilayer ribbon gaps have considerable uncertainty. We 
believe that the study is nevertheless useful because it 
can provide a framework for assessing the significance of 
future experimental findings. 

Our paper is organized as follows. We first summarize 
the DFT calculations that we have performed, comment- 
ing on the motivation for choosing various different semi- 



local approximations for the exchange-correlation poten- 
tial in section II. Then, in section III, we summarize the 
results that we have obtained for bilayer ribbon gaps, fo- 
cusing most extensively on the interplay between ribbon 
width, ribbon edge magnetism, and the external electric 
field between ribbon layers. Finally we summarize our 
results and present our conclusions. 



II. DENSITY FUNCTIONAL THEORY 
CALCULATIONS 

Our electronic structure calculations were per- 
formed with plane wave basis sets and ultrasoft 
pseudopotentials22. As an initial test to reproduce bulk 
bilayer graphene electronic structure and equilibrium in- 
terlayer separation, we placed the bulk bilayer graphene 
in a supercell with 10 A vacuum regions inserted along 
the direction perpendicular to the ribbons to avoid in- 
tercell interactions. The atoms were then fully relaxed 
without constraints. A 21 x 21 x 1 k-point mesh in 
the full supercell Brillouin zone (FBZ) was used with 
a 30 Rydberg kinetic energy cut-off. We estimated that 
these values yield a total energy converged to within 0.01 
meV/supercell. When combined with a LDA, these cal- 
culations yielded an interlayer separation of ~ 3.34A, 
within 1% of the experimental interlayer separation (3.35 
A). The same calculations with PBE96 and PW91 po- 
tentials overestimated the equilibrium layer separation 
rather badly (4.43 A and 3.65 A respectively). It is 
known that LDA or GGA does not include the van der 
Waals dispersion interactions as a result interlayer bind- 
ing energy and spacings are not expected to match with 
experiments 30 . The LDA success in predicting the in- 
terlayer distance is, therefore, not completely suprising 
as it has been shown previously 3 -^ that some fortuitous 
cancellation errors may be responsible for its success. 

For ribbon calculations, we introduced an additional 
transverse vacuum region of 10 A. We used 9x3x5 
k-point mesh in the FBZ and 30 Rydberg kinetic energy 
cut-offs in all cases. Total energy convergence was tested 
by using a larger k-point mesh, larger vaccum regions, 
and larger energy cut-off values. The electric field was 
applied perpendicular to the bilayer ribbons. We chose 
to apply several values of the electric field, including the 
value that was used in a recent experiment^ which was 
about 0.17 V/nm, and up to the maximum value close to 
the Si02 dielectric breakdown field of 1 V/nm 32 . Ribbon 
widths as large as 5 nm were considered. The c-orbitals 
along the ribbon edges were saturated with atomic hy- 
drogens. 

We do not consider edge roughness in our calculations, 
although there is a hint of roughness at the atomic scale 3 - 3 - 
in epitaxial graphene (grown on SiC substrate) and the- 
oretically, it was shown that considering edge roughness 
can have considerable effect on the electronic transport 
in single layer armchair ribbons 3 -^. We note that alter- 
nate edge functionalizations of the bilayers, other than 
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by hydrogen atoms, can also change the electronic struc- 
ture of bilayer zigzag ribbons? - because the localized edge 
states can react with different radicals thus altering the 
electronic structure of the ribbons. Our predictions are 
therefore most relevant to nanoribbons cut in a hydrogen 
environment and without significant edge disorder. 

Because we find that the LDA interlayer separation 
is close to the experimental value, we used only the 
LDA for our armchair ribbon electronic structure cal- 
culations. We relax the carbon atoms and the C-H dis- 
tances (initially chosen equal to the C-H bondlength in 
CH4 molecule) with the force threshold of 0.1 meV/A. 
The relaxed interlayer distances in armchair ribbons were 
found to be very close to the bulk values. As in the 
single-layer cas o 8 ' 36 , we identified three classes of arm- 
chair ribbons for both types of edge alignments. All bi- 
layer armchair ribbons were found to be non-magnetic, 
independent of the exchange-correlation approximation. 
(We refer to PW91 potential as the GGA below.) 

Zigzag ribbons with a-alignment led to non-magnetic 
ground states (with the initial configurations either as 
interlayer ferromagnetic or antiferromagnetic) when the 
LDA was employed whereas a magnetic ground state, 
which differs from the non-magnetic ground state in to- 
tal energy by only a few tens of meV, was obtained with 
the GGA. For /3-aligned ribbons, both LDA and GGA 
predict a magnetic ground state although GGA shows a 
stronger tendency toward magnetism. This means within 
DFT, the occurrence or absence of magnetism in bilayer 
ribbons is sensitive to the choice of the semilocal approx- 
imation. We stress here the fact that, in both a- and 
/3-aligned ribbons, the total energy difference between in- 
terlayer antiferromagnetic and ferromagnetic order is not 
large enough to call for a distinct magnetic state. Since 
our focus was to explore the broken symmetry states in 
zigzag ribbons and there was considerable uncertainty in 
predicting the magnetic order with GGA, for consistency 
and comparisons, we chose interlayer ferromagnetic order 
in both edge alignments and GGA as a semilocal approx- 
imation. Since with the relaxation of atoms the bulk in- 
terlayer distance was overestimated with GGA compared 
to the experiment, we chose to fix the interlayer separa- 
tion and the atomic coordinates at the experimental value 
for zigzag ribbon calculations. 

We note that the ferromagnetism which occurs in 
single-layer zigzag ribbons is thought to be related^ to the 
flat band which appears in TB models and splits when 
edge magnetic order is allowed. It is not surprising that 
adding the second layer can destroy the magnetic order 
since it also splits the flat band. As we show later in this 
article, a flat band also appear in a-aligned ribbons but 
is shifted from the Fermi level whereas the flat band lies 
at the Fermi level in /3-aligned ribbons. The position of 
the flat band with respect to the Fermi level may explain 
the magnetism in one alignment but not in the other. 



III. GRAPHENE BILAYER RIBBON GAPS 

We now present our results for the width, magnetism, 
and external electric field dependence of the bilayer gaps 
in ribbons with both edge type and the edge alignment. 
Wherever possible we will compare the GGA zigzag rib- 
bon results with corresponding LDA results. 
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FIG. 2: (Color online) Variation of energy gap with ribbon 
for balanced bilayer armchair ribbons, (a) a- alignment and 
(b) j3- alignment. Three classes of ribbons, denoted as 3p, 
3p+l and 3p+2 are shown. For a-aligned ribbons, p = 1-13 
whereas for /3-aligned ribbons, p=l-7, 9, 13 



A. Balanced Armchair Bilayers 

In this section, we discuss the width dependence of 
gaps in armchair ribbons with both edge alignments. 
Figure [2ja) shows variations of the gap versus the rib- 
bon width for a-aligned armchair ribbons. Here we la- 
bel the classes by N = 3p, 3p+l and 3p+2 where p 
range from 1 to 13 which translates to ribbons with 
widths of 1-5 nm. As expected on the basis of previ- 
ous work, the ribbon- width dependence is smooth within 
three classes which are distinguished by the number N of 
carbon chains across the ribbon mod 3, with two classes 
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showing semiconducting and class showing a tendency 
towards metallic behavior. The metallic behavior does 
not appear in the corresponding single layer graphene 
DFT calculation^. We believe that the crossings of 3p 
and 3p+l curves may be due to the inability of DFT 
(with LDA or GGA) to predict the gaps accurately espe- 
cially for narrow gap ribbons. The metallic behavior in 
3p+2 ribbons may be ascribed, similarly, to DFT-LDA 
not being able to resolve extremely small gaps. As in the 
monolayer ribbons^, the gaps decrease with increasing 
width. 

For comparison and for illustrative purpose, we chose 
p = 3-7, 9, 13 for the /3-aligned ribbons. DFT again pre- 
dicts three classes of gaps (Fig. H£b)), with no metallic 
regime for wider 3p+2 ribbons. The gaps are found to be 
consistently larger compared to the gaps in a-aligned rib- 
bons. To give a typical example, the gap in a /3-aligned 
ribbon with the width of 3.32 nm (0.213 eV) is 50 % larger 
than the gap in the a-aligned ribbon with the width of 
3.20 nm (0.141 eV), although both contain N= 27 chains. 
Eventually, for very wide ribbons the bulk limit of a zero- 
gap semiconductor is approached. In comparison, the 
corresponding DFT energy gap for a monolayer ribbon 
with the width of 3.2 nm is about 0.3 eV— . The energy 
gaps of bilayer nanoribbons are in general smaller than 
those of monolayer nanoribbons due to the interlayer cou- 
pling. 

To understand qualitatively the origin of three classes 
of bilayer armchair nanoribbons in a-alignment, we ana- 
lytically solved for the energy eigenstate by transforming 
the Hamiltonian to the one-dimensional TB model of a 
coupled two-leg ladder system (Figj3j). We calculate the 
bands only at the ribbon T-point - at which all gaps are 
minimized. 

We made two assumptions: (1) only nearest neigh- 
bor (NN) intralayer hopping t and (2) the NN interlayer 
hopping tj_ is allowed. For /^-alignment this analysis can 
not be performed due to the dangling bond present with 
this particular edge alignment, but we believe that /3- 
alignment qualitatively follows the same tendency as the 
a-alignment. 
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FIG. 3: (Color online) Schematic illustration of the TB 
Hamiltonian for the coupled two-leg ladder system. The la- 
bels used here are explained in the text. 



In Figure [3j a solid line represents NN intralayer hop- 
ping t. We consider only NN interlayer hopping tj_, 
which links bottom layer bi sites to top layer b\ sites for 
i = 1, 2, • • • , N. Then the TB model gives 



ea n = i(6 n _i + b n+l ) + tb n , 

ea' n = tib'^ + b' n+1 ) + tb' n , 

eb n = t(a n -i + a n+ i) + ta n + t±b' n , 

eb 'n = *( a n-l + a 'n+l) + ta n + t l- b n- 



(1) 



Let us define 



v± — 



(a„ ± a' n )/s/2 and 0± = (b n ± 



b' n )/V2. Then Eq.© can be rewritten as 



t{Pn-X+Pn+l)+tPt, (2) 

Assuming a± ~ A ± e m0 and /3± - B ± e mf> , we get 



£ 

-2t cos( 



-2tcos6> - t 
t £Tt± 



0. (3) 



Thus energy spectrum is given by 



-±,± 



±(t±/2) ± v / (^/2) 2 + (2tcos0 r + <) 2 (4) 



where the edge boundary condition results 9 r = m/ (N + 
1) for r = 1, 2, • • • , N. Note that when N = 3p + 2, 
there are zero-energy states e~- + = = at r = 

2(p + 1). This means that as in the monolayer case, 
bilayer armchair graphene nanoribbons are metallic for 
N = 3p + 2, whereas for N — 3p and N = 3p + 1, they 
are semiconducting. 



B. Unbalanced Armchair Bilayer Ribbons 

For the three classes of the armchair ribbons, in both 
edge alignments, we now discuss the external field effect 
on the gaps. Interestingly, we find that for a-aligned rib- 
bons with gaps below ~0.2 eV, the electric field has the 
effect of increasing the gap whereas for those above ~0.2 
eV, the gap decreases with electric field, as shown in Fig- 
ure |U This can be understood by using second order 
non-degenerate perturbation theory to semiconducting 
ribbons and first order degenerate perturbation theory 
to metallic ribbons. Below we show that there exist a 
critical gap, given by e c gap — [(y/5 — l))/2]t±, that can 
explain the electric field effects. With t± = 0.34 eV— , 
we get e c gap = 0.21 eV. 

To understand the influence of an external electric po- 
tential on the gaps, we consider the external potential 
difference between the layers denoted as U as a small 
perturbation to the on-site energies. Then Eq.fl]) and 
Eq.([2]) are modified as 



(e-U/2)a n = i(6„_i + b n+1 ) + tb n , 

(e + U/2)a' n = t{b' n _ 1 + b' n+l )+tb' n , 

(e-U/2)b n = i(a„_i + a„+i) + ta n + t±b' n 

{e + U/2)b' n = t{fl' n _ i + a' 1 ^i)+ta' n + t ± b n 



(5) 
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FIG. 4: (Color online) Variations of energy gap of the three 
classes of bilayer armchair ribbons with applied electric field 
and widths between 1-5 nm for a-aligned ribbons. For clarity, 
the ribbon with the width of 4.92 nm is denoted by stars 



change of energy levels can be obtained by second-order 
perturbation theory as follows: 



5e\ 
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(8) 



If these energy levels are low-energy states near the Fermi 
energy the unperturbed energy gap is e° gap — 2(7-/3) 
and the energy gap change due to the perturbation is 
given by 



Sr — 



7 2 



Alternatively, we can rewrite Eq.@ using £® ap , 
t± as following form: 



SSgap — C(s 



gap 



-gap 



)U 2 



(9) 



U and 



(10) 



where C is a positive constant and e gap = [(\/5— 1)/2] t±. 

U decreases the energy gap, 



> £ 



gap- 



Note that if £° _ 

while for s^ ap < s gap , U increases the energy gap in sec- 
ond order. Thus an external potential can increase or 
decrease the energy gap depending upon the size of the 
gap in a-aligned bilayer armchair graphene nanoribbons. 

For metallic armchair nanoribbons, a = 2t cos9 + t = 
in Eq.© for low-energy states near the Fermi energy 
thus zero-energy states are degenerate. The Hamiltonian 
in the zero-energy subspace becomes 



II — I U/2 \ 
~ 1 U/2 J 



(11) 



and 



± U T 



± 
n ' 



(6) 



Assuming at ~ A ± e in9 and /3± ~ B ± e md , we finally 
obtain the following Hamiltonian problem: 



thus a small perturbation U opens an energy gap linearly 
as 5e gap = U. 

In summary, we find that for ribbons with gap below 
£ gap (Fig- Sib), the gap increases with the electric field 

it 



gup- 



whereas for those ribbons with the gap above e' 
decreases with electric field (Fig. 0Ja)). We note that 
the /3-aligned ribbon gap show similar tendencies under 
interlayer external electric fields (figure not shown). 



/ 2tcos6 + t U/2 \ 

2tcos8 + t t± U/2 I 

U/2 2tcos6 + t 

\ U/2 2tcos6 + t -t± J 

(7) 

in the (A + ,B + ,A ,B ) basis. W can achieve a quali- 
tative understanding of the dependence of bilayer ribbon 
gaps on an external potential by treating U as a pertur- 
bation. 

For simplicity, let's define a — 2tcos9 + t, 2/3 = t±, 
7 = \J a 1 + /3 2 and 5 = U/2. The unperturbed energy 
levels from the lowest value are given by = —7 — /3 
, el = -7 + , eg = 7 - /3 and e\ = 7 + (3. The 



C. Balanced Zigzag Bilayer Ribbons 

In this section we address bilayer zigzag ribbons, for 
which the physics of gaps cannot be separated from the 
physics of edge magnetism. With the LDA a-aligned rib- 
bons were found to be non-magnetic, using either ferro- 
magnetic and antifcrromagnetic arrangements as a start- 
ing configuration. On the other hand non-zero mag- 
netic energy (tens of meV but slightly more than KbT), 
was obtained using GGA. The character of the interlayer 
magnetic coupling in these solutions (ferro or antiferro) 
could not be determined since the total energies were too 
close (within KbT). 
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In our view this signals a situation in which the compe- 
tition between non-magnetic and magnetic states is too 
close for reliable DFT predictions. For comparison with 
the /3-aligned ribbons below, we chose solutions with fer- 
romagnetic order between the layers and antiferromag- 
netic order across the layer in discussing the width and 
external electric field effects on bilayer gaps. 




FIG. 5: Band structure of bilayer zigzag ribbons in q- 
alignment in which the edge atoms are (a) without any mag- 
netic order and (b) with ferromagnetic order between the lay- 
ers. The width of the ribbon is 1 nm and the Fermi level is 
set at zero. Clearly the energy spectrum is gapless in (a) but 
has a finite gap in (b) . 

For a typical width of I nm, we plot the energy spec- 
trum of the a-aligned ribbon (Fig. EJa)). A gapless 
structure with a flat band shifted away from the Fermi 
level is clearly seen and it occurs roughly at one third 
the distance from the edge of the BZ. When we allow 
for magnetic edges in the calculation, due to a different 
spin order on A and B sublattices along the edges, a gap 
appears at the Fermi level (Fig. [5fb)) . 




magnetic order along the edge atoms lifts the degeneracy 
and a gap opens up in the energy spectrum (Fig. E{b)). 

The magnetic versus non-magnetic scenario in a- and 
/3-aligned ribbons respectively can be understood by rec- 
ognizing the position of the dispersionless state with re- 
spect to the Fermi level in the energy spectrum (Figs. 
[5]and[S|). In /3-aligned ribbons, the occurrence of a flat 
band at the Fermi level in the non-magnetic ground state 
results in a large density of states. The system is, there- 
fore, unstable towards developing long-range magnetic 
order (Stoner's criteria for itinerant magnetism). When 
exchange interaction is allowed between the edge carbon 
atoms, the system gains energy by opening up a gap in 
the energy spectrum and by simultaneously reducing the 
band energy, thus favoring a magnetic ground state. In 
a-aligned ribbons, in the non-magnetic ground state, the 
localized state at the Fermi level is absent and therefore 
the system does not show any clear tendency towards 
magnetism even when the exchange interactions are in- 
cluded. 
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FIG. 7: (Color online) Variation of the energy gap with the 
width of bilayer zigzag ribbons in both edge alignments using 
both LDA and GGA. For a-alignments, LDA predicts a non- 
magnetic ground state and thereofore no energy gap in the 
energy spectrum. 



FIG. 6: Same as Figure 5 but for /3-aligned ribbons. 

In /3-aligned ribbons, magnetic ground states were re- 
alized with LDA but the magnetic energy was not large 
enough to distinguish reliably between non-magnetic and 
magnetic states. However, GGA predicts a magnetic 
ground state with a magnetic energy of a few hundreds 
of meV, thus we can say reliably that it is magnetic. 
Again the precise form of the magnetic order was found 
to be uncertain because of the small energy differences 
between ferromagnetic and antiferromagnetic alignment 
between the layers. The energy band structure for the 
/3-aligned non-magnetic ribbons is shown in Figure [SJa). 
We observe a flat band at the Fermi level which is also 
seen in single layer graphene calculations. Allowance for 



Interestingly, we also found that the magnetic energy is 
almost independent of ribbon width for both edge align- 
ments. This can be understood by analyzing the mag- 
netic moments of the edge atoms for each ribbon width. 
We found nearly the same local moments (~ 0.3 [1b for 
a-alignment and ~ 0.32 fig for /3-alignment) on the edge 
atoms for all widths, and this translates to nearly the 
same magnetic energy per edge atom (0.2 meV for a- 
alignment and 0.25 meV for /3-alignment) for all ribbon 
widths. Since the magnitude of moments on the edge 
atom is independent of the width of the ribbons, the de- 
crease of the gap with ribbon width, shown in Figure 
UJ can be attributed to weakening quantum confinement 
effects for the 7r-orbitals. 

We note here that for /3-alignment the magnetism is 
weaker in LDA compared to GGA. Also shown in the 
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same figure are the bilayer gaps obtained with the GGA. 
These gaps are consistently larger than those obtained 
for a-aligned ribbons. This can be understood by com- 
paring the projected magnetic moments on edge atoms 
of the a- and /3-alignments. We found that edge atoms in 
/3-aligned ribbons carry larger moments than the atoms 
of a-aligned ribbons. As a result, the gap is expected to 
be larger. In summary, the opening of the gap in the en- 
ergy spectrum of the magnetic ribbons is associated with 
having opposite magnetic potentials on opposite ribbon 
edges. Nevertheless trends in the gap with the width 
should be understood in terms of quantum confinement 
effect. 



D. Unbalanced Zigzag Bilayer Ribbons 
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FIG. 8: (Color online) Variation of the energy gap with ex- 
ternal electric field for zigzag bilayer ribbons with (a) a- 
alignment and (b) /3- alignment. 



We apply external electric fields perpendicular to the 
ribbons with both edge-alignments (up to the dielectric 
breakdown field of Si02). We find that the gap decreases 
with increasing electric field in both cases. Although not 
explicitly proved by us for zigzag ribbons here, we believe 
that the gap will decrease with increasing electric fields 
for zigzag ribbons with widths larger than the critical 
gap (Fig. Ha) and (b)). 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have studied the electronic proper- 
ties of armchair and zigzag bilayer graphene nanoribbons 
both with and without the external electric fields using a 
first principles DFT-based electronic structure method. 
This paper summarizes our results for bilayer ribbons 
with AB (Bernal) stacking with two different edge align- 
ments which we refer to as the a- and /3-alignments. We 
find three classes of armchair ribbons the origin of which 
we explained using an analytical TB calculation. We dis- 
cuss the variation of the energy gap with applied electric 
fields, on the basis of a perturbation theory. A critical 
value of the bulk energy gap exists which controls the 
sign of the observed behavior. Gaps increase (decrease) 
with electric field for a bulk energy gap below (above) the 
critical value. Magnetic order in zigzag bilayer ribbons 
is found to be sensitive to the details of the semilocal 
exchange-correlation approximation. The local moments 
on edge atoms in zigzag ribbons are found to be very 
weakly dependent of the width of the ribbon, which im- 
plies that the gap dependence on ribbon width be purely 
consequence of quantum confinement of Tr-orbitals. By 
invoking band structure effects, we explained the mag- 
netic nature of the these ribbons. Although gap values 
are too sensitive to details of the exchange-correlation 
potential to allow fully predictive DFT results, we ex- 
pect that the present results will prove helpful in mov- 
ing toward a full understanding of how the experimental 
properties emerge from the interplay of magnetic order, 
width and the external electric field strength. 
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